Electrostatic confinement of electrons in graphene nano-ribbons 
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Coulomb blockade is observed in a graphene nanoribbon device with a top gate. When two in- 
junctions are formed via the back gate and the local top gate, electrons are confined between the 
pn-j unctions which act as the barriers. When no pn-j unctions are induced by the gate voltages, 
electrons are still confined, as a result of strong disorder, but in a larger area. Measurements on five 
other devices with different dimensions yield consistent results. 
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Confinement of the Dirac particles is of particular im- 
portance for the realization of nano-electronic devices in 
graphene such as quantum dots [1]. These would en- 
able one to perform single-level spectroscopy of Dirac 
particles, study their spin and valley degrees of freedom, 
and explore their potential for quantum coherent con- 
trol In conventional semiconductors, particles can be 
confined by potential barriers created via electrostatic 
gates. This approach permits independent control of the 
number of electrons on the island, the tunnel coupling be- 
tween the island and the reservoirs, as well as the tunnel 
coupling between neighbouring islands. Such flexibility 
and versatility has been instrumental for a wide variety 
of mesoscopics experiments. In graphene, this approach 
normally fails, due to the absence of a bandgap and the 
presence of Klein tunnelling [3, 4]. In previous studies, 
graphene has been etched into small islands, separated 
from the reservoir by narrow constrictions [5|, |6|, [fl] , but 
here it is difficult to tune the barriers. Alternatively, a 
bandgap could be created in graphene first, so that elec- 
trostatic gates can again be used for confinement. Theo- 
retically, a bandgap is predicted in graphene nanoribbons 
(GNRs) due mainly to quantum confinement [sl, Isl, lid ]. 
Experimentally, a transport gap has indeed been ob- 
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served in GNR devices 
is still under debate. 

Here we experimentally investigate GNR devices with 
a local top gate (TG) and a global back gate (BG) where 
the transport gap in the GNR enables electrostatic con- 
finement by the gates. Electrons are confined in an is- 
land where the barriers are formed by the pn-j unctions 
induced at the two edges of the TG, as demonstrated 
by the capacitances analysis of the measured Coulomb 
blockade. On the other hand, when no pn-j unctions were 
induced by the gates. Coulomb blockade was also ob- 
served, showing a larger confinement area. Here the is- 
land may be due to Anderson localization. Consistent 
results were found in five other devices with different di- 
mensions. 

Six devices {A to F) are fabricated on graphene flakes 
deposited on a substrate by mechanical exfoliation of nat- 
ural graphite [lH. The substrate consists of highly p- 




FIG. 1: (Color online) (a) Scanning electron microscope im- 
age of device E (scale bar 300 nm). The green dashed lines 
indicate segment I and II discussed in the text, (b) Schematic 
drawing of the device, (c) Current as a function of the back 
gate and top gate voltages for device A at Vhias = 150 /xV and 
T = 4 K. The insets illustrate the potential landscapes that 
are created in the four corners of the plot. They represent en- 
ergy diagrams along the ribbon length, where the gray band 
indicates the transport gap and the dashed lines represent the 
Fermi level Ep- A^/(//) indicates the doping level in segment 

I (II). 



doped silicon, acting as a back gate (BG), capped by 
285 nm of Si02. From their optical contrast against 
the substrate, we estimate that the flakes are single-layer 
[lit • Three electron beam lithography steps were used for 
patterning the devices. First, selected graphene flakes 
are patterned into GNRs, using PMMA as an etching 
mask and an Ar plasma for etching (for device F an O2 
plasma was used). Next we pattern a single top gate 
across each ribbon. The TG consists of 10/5/40 nm 
thick evaporated Si02/Ti/Au, and it covers only part 
of the ribbon, denoted as segment I (see Fig. la). The 
remainder of the ribbon (segments II) connects to wider 
pieces of graphene, which are contacted by 10/40 nm 
thick Ti/Au source (S) and drain (D) electrodes. The 
device is schematically illustrated in Fig. lb, and the 
relevant device dimensions are given in Table [D 

All measurements were performed in a ^He system at 
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GNR width 
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60 


50 


50 


50 


40 


40 


GNR length 
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2000 


1500 


1000 


700 


520 


520 


TG width 
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500 


400 


200 


140 


100 


50 



TABLE L The dimensions of devices Ato F. The TG dielec- 
tric is d = 10 nm thick Si02 for all devices. 



a base temperature of 350 mK, unless stated otherwise. 
We measured the two terminal resistance through the top 
gated GNR devices by applying a DC voltage bias, Vbias^ 
on the source electrode and measuring the current at the 
drain electrode. 

By tuning the BG and TG voltages, we can shape the 
potential landscape along the ribbon. Fig. Ic shows the 
low bias conductance of device A as a function of Vbg 
and Vtg at T = 4 K. Along the dark vertical band, the 
conductance is suppressed as Ep is within the transport 
gap in segment IL Along the dark diagonal band, the 
TG and BG dope the graphene with opposite polarity 
and Ef lies in the transport gap in segment L At zero 
gate voltages, the device was unintentionally hole doped. 

In the lower right (upper left) corner of Fig. Ic, the 
ribbon is in a npn {pnp) configuration. In this regime, 
holes (electrons) can be confined in the area (segment I) 
between the two pn-junctions owing to the presence of the 
transport gap [l7|, ll8| . We thus expect Coulomb blockade 
in the npn and pnp regimes. In the lower left (upper 
right) corner, the ribbon is in d^pp'p {nn'n) configuration. 
Here Fabry-Perot type resonances could occur between 
the two steps in the potential landscape, but no Coulomb 
blockade is expected in an ideal ribbon, as there are no 
barriers. The difference in energy from Ef to the middle 
of the transport gap in segment I and II are denoted as 
AEj and AEu^ respectively, which is a measure of the 
doping level. 

In the npn configuration, we observe pronounced cur- 
rent peaks separated by zero-current regions as Vtg is 
swept. A representative measurement is shown in Fig. 
2d for device 5, measured in the gate voltage configu- 
ration indicated by the black arrows in Fig. 2a and 3b. 
High bias measurements in the same range (Fig. 2e) 
show diamond-shaped regions in the Vbias — Vtg plane, 
in which current is blocked. Both are characteristic of 
Coulomb blockade due to the formation of an island that 
is only weakly coupled to the leads. 

In this device, over 700 Coulomb peaks were resolved in 
the range —2 V< Vtg < —0.2 V, corresponding to a large 
change in doping level in segment I, —360 meV^ AEj ^ 
—240 meV (Fig. 3a). Here AEj reflects the doping level 
in segment I (Fig. Ic inset), and is roughly estimated 
by considering the density of states of bulk graphene, 
AEj = {±)TiVF^/^ni^ where n/ is the carrier density 
in region I, vf is the Fermi velocity of bulk graphene, 
and the (+) and (— ) signs represent electron and hole 



doping, respectively. The spacings between neighbouring 
peaks, AVr^Q^ ^ are shown in Fig. 3a (black spheres), 
as a function of the peak positions. The average value 
(AVrPo^) = 2.0±0.4 mV corresponds to a TG capacitance 
Ctg = 70 — 100 aF, close to what one would expect 
from simple parallel plate capacitance between the TG 
and segment I, C^^q = coCrWl/d = 70 aF, where = 
3.9 is the relative permittivity of Si02. In addition, the 
capacitance to the back gate, measured to be ~ 3.9 aF, 
compares well to the value expected from the geometry 
of an island of area {wl). The agreement demonstrates 
that for this device, an island is formed between the two 
pn-j unctions in the npn configuration. 

In addition, we measured over 100 Coulomb diamonds 
similar to Fig. 2e, and the extracted addition energy 
■^npn each diamond is shown in Fig. 3a with green 
triangles. The average addition energy is (E^^^) = 1.0± 
0.4 meV. From Fig. 3a, no shell filling or evident top gate 
voltage dependence is observed in either AVtg or but 
both quantities show a large spread similar to [1, Q , due 
to contributions from both the level spacing and strong 
disorder from the ribbon edges, which is discussed further 
below. In all these measurements, segments II of the 
GNR were heavily n-doped such that disorder in the leads 
was largely screened {Vbg = +81 V) 0. 

Unexpectedly, Coulomb blockade was also observed 
when no pn-j unctions are present. Fig. 2b and 2c 
show a representative current trace and Coulomb dia- 
monds measured from the same device in a pp^p config- 
uration indicated by the red arrows in Fig. 2a and 3d. 
Over 1400 current peaks were observed in the range —1 
V< Vtg < 2.2 V, corresponding to a change in dop- 
ing level AEi from -430 meV to -230 meV (Fig. 3c). 
We measured 70 Coulomb diamonds in the same range. 
The extracted peak spacings AVr^Q^ and addition energy 
E^'P P are shown in Fig. 3c with red spheres and blue tri- 
angles, respectively. The measurements were taken at 
Vbg = (red arrow in Fig. 2a), such that segment II of 
the GNR was heavily p-doped. 

The average peak spacing in the pp'p configuration 
{AV^^P) = 1.8 ± 0.4 mV is very close to (AF^g^), but 
the average addition energy {E^^ ^) = 0.5 ± 0.2 meV is 
only half the value of {E'^^^). The back gate capaci- 
tance in pp^p is ~ 10 aF, indicating an island area of 
50 nm by 700 nm, larger than {wl) (assuming the island 
extends over the entire ribbon width in the transverse 
direction). All average quantities were reproducible over 
multiple thermal cycles. Therefore the island formed in 
the pp'p configuration is located in part under the TG, 
but extends to a larger area than the island in the npn 
case. 

Coulomb blockade in the npn configuration was found 
in three other devices of different dimensions (A, C, E) 
where the npn regime could be accessed (the various de- 
vices exhibited different positions of the charge neutrality 
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FIG. 2: (Color online) (a) Current as a function of Vbg for 
device B at Vtg = -1.2 V and Vbias = 200 /xV. The top 
axis indicates the corresponding doping level in segment II, 
AEii, estimated in the same way as AEj (see text), (b) 
Coulomb oscillations as a function of Vtg in the pp'p config- 
uration. Vbg = as indicated by the red (gray) arrow in (a) 
and Vbias = 100 /xV. (c) Differential conductance dl/dV as 
a function of Vtg and Vbias (Coulomb diamonds), measured 
in the same ppp regime as (b). (d) Coulomb oscillations as 
a function of Vtg in the npn configuration. Vbg = +81 V 
as indicated by the black arrow in (a), and Vbias = 100 /xV. 
(e) Coulomb diamonds measured in the same npn regime as 
(d). The addition energy Ea is taken as the average of the 
two blue arrows. 



point). For device also the pnp configuration could be 
reached, where the measurement results are analogous 
to those for npn. The extracted (AF^g'') and {E^p"") 
for these devices are summarized in Fig. 4a (black filled 
circles) and 4b (black filled squares), respectively. For 
all four devices, the measured {/SV^q^) agrees quantita- 
tively well with the numerically computed c/Ctg by 
considering an island of size vol (black dotted line in Fig. 
4a). Moreover, both {AV^q^) and {E^^^) increase with 
decreasing area iwl)^ consistent with the observation in 
device 5, i.e. carriers are confined in segment I in the 
npn configuration. 

Furthermore, we observed Coulomb blockade in the 
pp'p regime in all six devices Aio F. For devices A^B^E, 
we could also access the nn'n configuration and the re- 
sults are similar to those in the pp^p case. For all devices, 
the average addition energy, {E^^ ^) in the pp'p regime 
is much smaller than that in the npn case, and does not 
vary much despite the differences in device dimensions 
(except for device F), as shown by the red open squares 
in Fig. 4b (no clear Coulomb diamonds were observed in 
device A in the pp'p regime). In devices B and the 
measured back gate capacitance indicates that the island 
sizes are 50 x 700 nm^ and 50 x 250 nm^, respectively 
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FIG. 3: (Color online) (a) A Vtg (black spheres) and Ea 
(green triangles) as a function of peak positions in the npn 
regime for device B. Vbg — +81 V. (b) Current (black line) 
as a function of Vtg at Vbias — 100 jiY and Vbg — +81 V. 
(c) AVtg (red spheres) and Ea (blue triangles) in the pp'p 
regime for at Vbg =0. (d) Current (red Une) as a function 
of Vtg at Vbias = 100 /xV and Vbg = 0. The corresponding 
doping level in segment I, A^j, are indicated in the top axis 
(a) and bottom axis (d). The two insets illustrate the po- 
tential landscapes created at the gate voltage configurations 
where Fig. 2b and 2d are taken. 



(Fig. 4b inset), significantly larger than the area of seg- 
ment I. The capacitance of this large island to the relative 
narrow TG, Cj^'^, is still roughly the same as that for 
an island limited to segment I, so the peak spacings in 
pp'p (Fig. 4a, red circles) are similar to those for npn. 

The reproducible scaling of addition energy and peak 
spacing as a function of devices dimensions is consis- 
tent with the results obtained from device B: when 
pn-j unctions are induced by the TG and BG, an is- 
land is formed in between the junctions; without the in- 
junctions, a much larger island is formed, presumably due 
to disorder. However, the source and drain capacitances 
are comparable {npn) or even larger {pp'p) than Ctg^ 
and contribute more than half of the total capacitance. 
This means that extracting the island size from Ea may 
be unreliable. For device B {npn d^ndpp'p) and D {pp'p), 
we have measured the back gate capacitance, which gives 
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an independent estimate of the island size (Fig. 4b inset) 
and is in agreement with our interpretation of the island 
size in the two regimes. 

If islands are induced by disorder in the pp'p case, there 
are likely to be disorder-induced islands in the npn case as 
well in addition to the islands formed by the pn-j unctions. 
Indeed evidence of multiple islands was observed exper- 
imentally in several cases (not shown). The presence of 
the additional islands contributes to the large spread in 
peak spacings and addition energies mentioned earlier. 

These disorder-induced islands in the pp'p regimes are 
5-10 times longer than the ribbon width, which could 
be explained by Anderson localization, due to strong 
scattering at the rough ribbon edges as proposed in 
0, El, 113, This reveals a different aspect of 

the electronic properties of the GNRs compared to other 
work, where the extent of the island is found to be compa- 
rable with the ribbon width [14, 2^, [l^]. Further studies 
are needed in order to clarify the underlying mechanisms 
behind the various observations. 

In conclusion, a single electron transistor is formed in 
graphene nanoribbon devices with single top gates. Two 
pn-j unctions at the two edges of the top gate induced by 
the top gate and back gate voltages act as barriers to 
form an island. Hundreds of Coulomb peaks were ob- 
served in this regime. In the absence of the pn-j unctions, 
regular Coulomb blockade is also observed where the is- 
land can be induced by ribbon edge disorder. Observa- 
tions from measurements of five other devices give consis- 
tent results. We anticipate that multiple top gates on a 
graphene nanoribbon will offer additional control for fu- 
ture device applications, and provide further insight into 
the electronic properties of graphene nanoribbons. 
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FIG. 4: (Color online) (a) Average peak spacing versus the 
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in the npn (black filled squares) and ppp (red open squares) 
configurations versus {wl)~'^. The error bars in (a,b) repre- 
sent standard deviation. Inset: the island size in the direction 
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The underlying energy diagrams for the nnn (c) and pnp (d) 
regimes are also shown where the gray band represents the 
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